Decoherence and Matter Wave Interferometry 
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A two-slit interference of a massive particle in the presence of environment induced decoherence 
is theoretically analyzed. The Markovian Master equation, derived from coupling the particle to 
a harmonic-oscillator heat bath, is used to obtain exact solutions which show the existence of an 
interference pattern. Interestingly, decoherence does not affect the pattern, but only leads to a 
reduction in the fringe visibility. 
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I. INTRODUCTION 



Recently it has been possible to observe diffraction of 
large molecules like Cqq, which are expected to behave 
like classical particles Interference effects have also 
been observed in larger molecules [|J . These effects are a 
consequence of the linear superposition principle applied 
to the wave functions of the particles and can only be de- 
scribed quantum mechanically. For massive particles, the 
existence of such superpositions is classically unimagin- 
able and is not a familiar observation in the real physical 
world . The quantum-classical transition and the nature 
of classicality as an emergent property of an underlying 
quantum system has been the subject of a lot of study 
in recent years 0, 0| ■ The decoherence phenomenon has 
been widely discussed and accepted as the mechanism re- 
sponsible for the emergence of familiar classical features 
in the real physical world. Decoherence results from the 
irreversible coupling of the system to its environment. 
The recent experimental observation of diffraction and 
interference patterns for large molecules raises a natural 
question: How far can one go before decoherence effects 
destroy the interference pattern of massive objects? The 
effect of decoherence on matter wave interferometry has 
been explored quite well experimentally. Hornberger et 
al have studied the effect of decoherence due to collisions 
of the interfering particle with gas molecules They 
observed that with increasing pressure of the gas, the 
fringe visibility goes down. Chapman et al have studied 
the effect of photons scattering off interfering atoms [6| . 
Kokorowski et al have studied the same effect using mul- 
tiple photons [?J • Sonnentag and Hasselbach have studied 
effect of decoherence in electron biprism interferometer 
Q. On the theoretical side there have been several fo- 
cussed studies on specific systems [9[[Toj. However, there 



is a need to understand the effect of decoherence on the 
observed interference pattern in a double-slit experiment 
in a simplified way. It is intuitively obvious that the 
resulting interference pattern in such a scenario would 
be a consequence of the interplay between the strength 
of decoherence, the slit separation and the distance the 
particle travels from the slit to the screen. One would like 
to understand to what extent each of these parameters 
play a role. 

Decoherence effects on quantum superpositions have 
generally been studied by assuming an initial state which 
is a superposition of two spatially localized wave-packets. 
For an initial state of the form ip(x) = ipi (x) + ip2 {%), the 
density matrix in the position representation is 



p{x,x') = Mx^Kx^ + Mx^Ux') 



(1) 
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For localized wavepackets (which best describe a massive 
particle), the four terms in this density matrix correspond 
to four peaks. The last two terms correspond to the off- 
diagonal peaks. Decoherence arguments show that the 
off-diagonal peaks of the density matrix die out, in time, 
because of the effect of the environment [4| . Thus, the ap- 
pearance of classical behaviour via decoherence is marked 
by the dynamical transition of the pure density matrix 
to a statistical mixture. For example, using the Marko- 
vian Master equations with some approximations, Zurck 
has argued that the density matrix for a free particle in 
an initial coherent superposition of two Gaussian wave 
packets separated by Aa; decoheres (i.e. the off-diagonal 
peaks decay) over a time scale which goes inversely as the 
square of the separation, Ax 2 , between the two parts of 
the superposition Q ■ For classical systems and standard 
macroscopic separations, this decoherence time is shown 
to be almost 10 -40 times smaller than the thermal relax- 
ation time of the system. Macroscopic superpositions, 
are, thus, almost instantaneously reduced to a statisti- 
cal mixture, a situation which is classically interpretable. 
The decoherence approach has been used to study many 
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models in the context of quantum measurement and the 
decoherence mechanism has been explored in the experi- 
mental regime also [TH [T3| . It is now generally accepted 
that the two main signatures of the decoherence mech- 
anism are, (i) in the classical regime decoherence takes 
place over a time scale that is much smaller than the 
thermal relaxation time of the system, and (ii) the deco- 
herence time goes inversely as the square of the separa- 
tion between the two parts of the superposition [J, [l3| . 
If one were to look at a two-slit interference situation, 
one can write a state of the form ip(x) = ipi(x) + ip2(x), 
where ip\ {x) and ip2 (x) correspond to the probability am- 
plitudes for the particle to pass through slit 1 and slit 2, 
respectively. As is well known, the interference pattern 
corresponds to the position probability distribution of the 
time evolved wave function: 

= p(x,x,t). (2) 

One might be naively tempted to look at the off-diagonal 
components of the density matrix at the screen. However, 
the interference pattern is obtained from the probability 
distribution of of the particle on the screen, which is just 
the diagonal components of the density matrix. So, it is 
not obvious if the dying out of off- diagonal components 
of the density matrix, at the screen, as a consequence 
of decoherence also corresponds to a disappearance of 
the interference pattern. There is a need, therefore, to 
study the evolution of the state of the particle along with 
the effect of the environment, and analyze the emerging 
position probability distribution for the existence of an 
interference pattern. Savage and Walls have addressed 
this issue by studying the evolution of a superposition of 
two plane-waves under the influence of decoherence, and 
its effect on the interference 14]. However, though their 
results illustrate the effect of decoherence on the inter- 
ference pattern, their use of plane waves to describe the 
state of the interfering molecules seems a little unrealis- 
tic as it is obvious that large molecules would be better 
described by localized wave-packets. 

Recently, decoherence effects on two-slit diffraction has 
also been theoretically analyzed by treating the effect of 
the environment using certain phenomenological models 
fT?| . However, a fully quantum mechanical analysis, 
using a microscopic model of the environment has not yet 
addressed this issue. This kind of analysis is very rele- 
vant from the point of view of quantitatively probing the 
elusive quantum-classical boundary. The study is also 
highly relevant in the light of several recent proposals to 
exploit purely quantum mechanical features for quantum 
computation and quantum information processing [l5|. 

In the following, we present an analysis of matter-wave 
interferometry in the presence of a dynamical, quantum 
environment. Starting with an initial superposition de- 
scribing the double slit situation, we study the dynamics 
of the system through the Markovian master equation 



which takes into account the coupling to the environ- 
ment. Exact solutions for the position probability distri- 
bution clearly bring out the role of environment-induced 
decoherence on the interference pattern. The rest of the 
paper is organized as follows. In Sec. II we set up the 
theoretical frame work for our analysis of decoherence 
and present our results. Further, we discuss fringe visi- 
bility and some of the factors affecting it. Finally, in Sec. 
Ill we summarize the main conclusions of this work. 



II. THEORETICAL ANALYSIS 

A. Coupling to the environment 

The effect of decoherence on the quantum evolution of 
a system can be studied by coupling it to a model envi- 
ronment. A popular model for the environment is a set 
of non-interacting harmonic oscillators, which may arise 
out of different physical situations. The Hamiltonian for 
a "free" particle, coupled to such a model environment 
can be written as 

k K 

Here x and p denote the position and momentum of the 
particle of mass to, and the second term represents the 
Hamiltonian for the bath of oscillators (environment). 
Xfc and Pfc are the position and momentum coordinates 
of the kth harmonic oscillator of the bath, C^s are the 
coupling strengths and f^'s are the frequencies of the os- 
cillators comprising the bath ,18]. As one is interested in 
the dynamics of the particle, and not the detailed dynam- 
ics of the environment (which is not within one's control 
anyway) , it is conventional to look at the reduced density 
matrix of the system, where environment variables have 
been traced over. 

This Hamiltonian has been studied extensively, to ob- 
tain the dynamics of the reduced density matrix. Of 
particular relevance here is the case of "ohmic coupling" , 
which gives the correct limit of classical dissipation. For 
our analysis, we deal directly with the Markovian master 
equation for the reduced density matrix of the system 
in the position representation, where the environmental 
degrees of freedom have been traced out[H, liH [20j: 

dp r (x,x') ih ( d 2 p r d 2 p r \ 
dt ~ 2ro [jkc 2 ~ ~dxJ 2 ) 

, ( dp r dp r \ D 2 

-^ X - X) \-dx--M)-^ { *-*)Pr- 

(4) 

Here D = 2m r ykBT for a thermal bath. Equation ([J| 
has been derived by assuming the harmonic oscillator 
environment to be in equilibrium, at a temperature T. 
The parameter 7 can be assumed to signify the strength 
of the coupling of the particle to the environment - it has 
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its origin in the coupling strengths Ck appearing in the 
Hamiltonian ([3]). 

At this stage it might be worthwhile to point out that 
this model of dccoherence is fairly general, and spans a 
wide range of physical situations. For example, dccoher- 
ence due to a photonic heat bath can be described by 
a particle interacting with the modes of the electromag- 
netic field, modeled by harmonic oscillators. In the case 
of intereference of electrons passing close to a conducting 
plate H|, the decoherence can be described by bosonic 
excitations of a Fermi sea of conduction electrons - this 
again, can be done easily by using a harmonic oscilla- 
tor heat bath. In addition to this, the Master equation 
(j4|) can also be arrived at by studying the the quantum 
evolution of a particle undergoing collisions with smaller 
particles, using a stochastic formalism [201 ] . This relates 
very well to the experimental study of collisional deco- 
herence in matter- wave interferometry [2lj . 



B. Decoherence 

Now that the framework of our analysis is set up, let 
us get back to matter wave interferometry. The particle 
encounters the double-slit, after which it travels a dis- 
tance L, say, along the y-direction, before registering on 
the screen. Clearly, for the interference pattern to be visi- 
ble, coherence along the x-direction is important, whereas 
the dynamics along the y-axis just serves to transport 
the particle from the slits to to the screen. In order 
to simplify the calculations, we assume that the parti- 
cle emerges from the double-slit, travels along the y-axis 
with a well-defined average momentum po and reaches 
the screen after a time tj, = mL/p Q . We now focus only 
on the time evolution of the particle's wave function in 
the x-direction. 

At this stage, let us specify the functional form of the 
initial state that emerges from the double-slit. We as- 
sume that the action of the slit is to prepare a super- 
position of two Gaussian wave-packets, centered at the 
location of the respective slits, with a width equal to the 
width of the slit. We define the initial state as 



#r,0) 



1 



1 



~(x-x ) 2 /e 2 



-(x+x ) 2 /e' 



(5) 



where xq — d/2, d being the separation between the two 
slits. (5) represents two Gaussians centered at x = ±xq. 

It may be noted that we have ignored the decoherence 
effects that could occur before the particle reaches the slit. 
From the point of view of an analytical calculation, it is 
difficult to take into account the effect of decoherence 
both before and after the double slit. Moreover, there 
could be specific experimental situations where decoher- 
ence occurs only after the double slit. In the experiment 
of Sonnentag and Hasselbach on interference of electrons, 
the decoherence takes place due to a metallic plate kept 
after the double-slit 



For the time evolution after the double slit, the initial 
density matrix for the particle state can be written using 
©, as 



p(R,r,0) = 



1 



V27 



(6) 



where R — x + x' and r — x — x' and d is the slit sep- 
aration. In order to analyze interference, one needs to 
obtain the density matrix of the particle at a time tL, 
when it reaches the screen. For this, we solve the master 
equation ([?]), with the initial condition It turns out 
that an exact solution for p r (x, x, ti,) can be obtained. 
This will contain full information about the interference 
and the degree of decoherence. However, for this purpose 
we do not need the full density matrix of the particle - 
we are only interested in the position probability distri- 
bution of the particle on the screen, which is given by 
the diagonal part [x = x') of the reduced density matrix. 
From the exact solutions one can see that the position 
probability distribution of the particle on the screen has 
the final form: 



p r (x,x,t L ) 



1 -(x+x ) : 



cos 



xdh(l 



+ 

27*1.) 



where T — 



16m 2 e 2 7 3 



{4~/t L + 4e~ 2 " /tL - e 



•47*1 



(7) 

3), n 2 = 



Starting from (|U) , the result ([7]) is exact and represents 
the position probability distribution at the screen. The 
presence of the cosine term indicates the existence of the 
interference pattern. However, (7) is a very general result 
which can be used to study many things, like damped 
motion of the particle, the effect of friction and of course, 
decoherence. In order to relate it to decoherence, this 
result should be analyzed in a suitable limit. This limit is 
set by the requirement that the interaction of the particle 
with the environment be so weak that the dissipative 
effects are not noticeable, only the decoherence is. This 
will correspond to the limit I/7 3> ti, i.e, when the 
relaxation time of the system is much much larger than 
the time scales over which the experiment is performed 
and the observations made. In the following we make 
this approximation: 



9/ 



2ht 



3m 2 



A^ 2 



2Dt\ 



DX 2 ,L 2 



3m 2 e 2 6tt 2 ^ : 



-t L - 



(8) 



Here, we have assumed that = mL/po, which implies 
that Ml/th = XfjL/2-K, Xj being the de Broglie wave- 
length. In the expression for f2 2 , we have made a further 
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approximation, e <C XdL/ire, i.e., the spread of the wave- 
packets, after travelling a distance L, is much larger than 
the original width. This will happen when the original 
width is very small. This is a realistic expectation since 
the slits should be narrow enough to let the two wave- 
packets spread and overlap with each other to lead to an 
interference pattern on the screen. 

Using the approximations J5]), © reduces to 



p r (x,X,t L ) 



1 



1 -1.-.Q) 



1 -(x + xq) 2 



7T(T \ 2 



exp(- 



tjj . f TTXd 

cos < . 

24t d ' \ \ d L 



(9) 



where a = X^L/ne and tb = 2m^^TcP ■ R- ecai l m g ear- 
lier results on decoherence, one can recognize td as the 
decoherence time of a superposition of two wave-packets, 
separated by a distance d, due to interaction with an en- 
vironment at temperature T, with a coupling strength, 
or relaxation rate 70, [l3|. 

Note that without the term exp(— ^ ), © represents 
the interference (position probability distribution) of two 
wave-packets of initial width e each. The expression also 
represents the interference pattern corresponding to a 
matter- wave of de-Broglie wavelength A<j, having trav- 
elled a distance L from the double slits. The decoher- 
ence effects come in only through the term exp(— 2 \t d )• 
This term affords a simple physical meaning - two wave- 
packets, separated by the slit distance d lose their coher- 
ence after a time 775, which, by definition, is the deco- 
herence time. If the wave-packets reach the screen at a 
time tx, which is much much larger than td, no interfer- 
ence will be visible. For the interference pattern to be 
observable, td should obviously be of the order of t^- 

Since our analysis is not tied to any one experimental 
situation, we use the following arbitrary values for various 
parameters: d = 1 /xm, e = 0.1 /im, = 5 x 10 _6 //m, L 
= 20 cm. We assume that the Gaussian width of a wave- 
packet emerging out of a rectangular slit of width w will 
be roughly w/2. The dimensionless parameter t^/Tp is 
a measure of decoherence. Using these values, expression 
© is plotted in FIG. 1. Note that the overall Gaussian 
profile is due to the finite width of the slits. If the slits 
were infinitely narrow, one would see a flat profile with 
narrow interference peaks. 



C. Fringe visibility 

One quantity of particular importance in matter-wave 
interferometry is the fringe visibility and the effect of 
decoherence on it. Conventionally, it is defined as 



V = 



In 



(10) 



where I rnax and I m i n represent the maximum and min- 
imum intensity in neighbouring fringes, respectively. In 
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FIG. 1: The interference pattern p r (x, x,tL) plotted against 
the position on the screen x, for e = 0.1/im, Ad = 5 x 10 _6 /xm, 
L — 20 cm, d — 1.0/im, and for various values of ti/ro. Solid 
line represents ti/ro = 4.0, dashed line represents ti^/ru = 
20 and the dotted line represents ti/rn = 60. 



reality, fringe visibility will depend on many things, in- 
cluding the width of the slits. For example, if the width 
of the slits is very large, the fringes may not be visible 
at all. If we focus on just the effect of decoherence on 
fringe visibility we can assume that the slits are so nar- 
row that we get an essentially flat background profile. 
This means that a will be so large that the functions 

-(x±x ) 2 

e ^ will have the same values at the points of max- 
imum and minimum intensity. Maxima and minima of 
will occur at points where the argument of the cosine 
is 1 and -1, respectively. Taking two such neighbouring 
points, the fringe visibility can be written as 



V 



exp( 



24t p ■ 



cosh( 



(11) 



where x n denotes the mean position of nth fringe. 
Clearly, the fringe visibility goes down as th becomes 
larger than td — 2m-yk B Td 1 • This can most easily hap- 
pen when either 7 or m becomes large. 

The expression for fringe visibility, (fTTj) . details its de- 
pendence on various physical parameters. It can be writ- 
ten in an expanded form as follows: 



V = 



cosh( 



2x n Xn 



(12) 



If the present microscopic model is compared with cer- 
tain stochastic models of the environment, the parameter 
7 turns out to be proportional to the rate of collision of 
the interfering particle with the smaller particles consti- 
tuting the environment 20]. If one tries to relate this 
model to collisional decoherence, as in the experiment of 
Hornberger et. al. the pressure of the gas in the 
chamber should be directly proportional to the collision 
rate, and therefore to the parameter 7 in our calculation. 
Eqn. (|12[) indicates an exponential decay of visibility 
with 7. This implies that the fringe visibility should go 
down exponentially with increasing pressure. This is in 
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broad agreement with the experimental findings of Horn- 
berger et. al. [a|, and does not depend on the fact that 
Hornberger et. al. use a Talbot-Lau interferometer in- 
stead of a double-slit. If one relates the calculation here 
to the experiment of Sonnentag and Hasselbach on inter- 
ference of electrons, the decoherence is due to the elec- 
tron's interaction with the metallic electrons of the plate. 
In this case, 7 is given by [1] 



7 



2 

e p 



32nmz 3 ' 



(13) 



where p is the resistivity of the place and z is the distance 
of the interfering beam from it. So, in this case too, 
the decoherence effects are expected to be stronger when 
electrons are closer to the plate, or z is smaller. This 
would mean, larger value of 7. 

Fringe visibility is also very sensitive to the separation 
between the two slits, d. It decays exponentially with 
the square of the slit separation. To get an idea of the 
magnitude of the effect, for example, if a particular slit 
separation gives a visibility of 60 percent, doubling the 
slit separation will reduce the visibility to about 13 per- 
cent. 

It is also clear from (fT2)) that fringe visibility goes down 
exponentially with the mass of the interfering particle. 
So, if one were to use a molecule, say, twice as heavy 
as C*6o, one would have to either reduce the temperature 
by the same factor, or decrease the pressure by the same 
factor, in order to get the same visibility as for Cqq. This 
is in tune with the general expectation that the interfer- 
ence arising from the quantum nature of particles in the 
double-slit scenario analyzed here will be more vulnera- 



ble to decoherence for more massive particles. In many 
situations, the coupling to the environment may also be 
related to the mass, and such simplified logic may not 
always be correct. 



III. CONCLUSION 

We have theoretically analyzed the effect of decoher- 
ence on a massive particle, in a two-slit interference set- 
up, using a quantum dynamical model of the environ- 
ment. Interestingly, the interference pattern does not get 
distorted as a result of decoherence. The effect of deco- 
herence is only to reduce the visibility of the interference 
fringes. The fringe visibility crucially depends on the de- 
coherence time of a superposition of two freely evolving 
wave-packets, initially localized at the centers of the two 
slits. Our results clearly demonstrate the two main sig- 
natures of the decoherence mechanism, namely (a) the 
decoherence time is much smaller than the thermal re- 
laxation time, and (b) the decoherence time is inversely 
proportional to the square of the "separation" between 
the two parts of the superposition. Our results also show 
that the fringe visibility goes down exponentially with the 
mass of the particle, with the temperature of the environ- 
ment, with the square of the slit separation and with the 
coupling strength of the particle to the environment. Fi- 
nally, the results are in good qualitative agreement with 
the matter-wave interferometry experiments carried out 
to study decoherence effects with large molecules and also 
with electrons. 
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